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We extrapolate the first moments of the generalized parton distributions using heavy baryon 
chiral perturbation theory. The calculation is based on the one loop level with the finite range 
regularization. The description of the lattice data is satisfactory and the extrapolated moments at 
physical pion mass are consistent with the results obtained with dimensional regularization, although 
the extrapolation in the momentum transfer to t = does show sensitivity to form factor effects 
which lie outside the realm of chiral perturbation theory. We discuss the significance of the results 
in the light of modern experiments as well as QCD inspired models. 



I. INTRODUCTION 



The study of hadron structure and, in particular, the spin structure of the proton is one of the most exciting 
challenges facing modern nuclear and particle physics. Through measurements of high energy scattering of electroweak 
probes one can measure the parton distribution functions (PDFs), which describe the number densities of partons 
with momentum fraction, x, in the nucleon. The study of the spin and flavor dependence of the PDFs has provided 
a wealth of data which has proven critical to our understanding of the structure of the nucleon [l| . 

In recent years, the extension of this concept to so-called generalized parton distributions (GPDs) has attracted 
enormous interest 0, 0] ■ The GPDs are related to the amplitude for deeply virtual Compton scattering where the 
initial and final photons have different momenta. It is widely anticipated that the GPDs should provide even more 
information concerning the internal structure of the nucleon. Of particular relevance to our present work is the link 
between the low moments of the GPDs and the angular momentum carried by each quark flavor within theproton Q ■ 
This is now widely considered essential to a satisfactory resolution of the famous proton spin crisis [H, [f| , because 
successful models suggest that a great deal of the proton spin is carried as quark orbital angular momentum @-H| . 

There have been many theoretical and experimental studies of the GPDs. On the theoretical side, much of the 
work has centered on the most effect ways to parameterize them llfjf. There has also been quite a bit of work on 
phenomenological models, such as the MIT and cloudy bag model |lH. [l2j], the constituent quark model [HI, flij . 
the NJL model [l5|, the light-front model [TH, 0], the color glass condensate model [l8[ and the Bethe-Salpeter 
approach [l9|, [2(|. On the experimental side, various groups have focussed on different kinematic ranges. ZEUS and 
HI measured the DVCS cross section for the first time, with in the very low range 1CT 4 < xb < 0.02 [2l|, [22j. 
COMPASS focussed on a little larger x B - from ~ 0.01 to ~ 0.1 [Hj]. The data from HERMES is in the range 
0.03 < xb < 0.35 0, H^. The study of the high xb domain, where valence quarks should dominate, requires high 
luminosity and a high energy electron beam. With the JLab 12 GeV upgrade, xb will range up to 0.7, while the 
current JLab data is in the range 0.15 < xb < 0.55 26-28|. 

The GPDs have a close relationship with the form factors. By integrating the GPDs with different powers of the 
momentum fraction x, the GPDs can be transformed into Mellin moments. There has been important work on the 
moments and form factors from various lattice QCD collaborations [29l - l3~i1 ] . as well as chiral perturbation theory 
[35l - [37l ] . As for many other observables computed in lattice QCD, current lattice simulations have been concentrated 
at quite larg e pion mass. While chiral perturbation theory is only expected to be convergent at quite low pion 
mass [33|, Hi|, in order to relate the simulations at large quark mass to experimental data, the lattice data has been 
extrapolated using covariant and heavy baryon chiral perturbation theory with dimensional regularization - e.g., in 
Ref. Ill- 
Based on the observation that all hadron properties show a slow, smooth variation with quark mass above ~ 
0.4GeV, suggesting that chiral corrections (pion loops) are strongly suppressed there [32] . an alternative regularization 
method, namely finite-range-regularization (FRR) has also been used to extrapolate the lattice data. It was first 
applied in the extrapolation of nucleon mass and magnetic moment [38:-40] . The remarkably improved convergence 
properties of the FRR expansion mean that lattice data at large pion mass can be described very well and the 
obtained nucleon mass at physical pion mass is close to the experimental value. Later, the FRR method was applied to 
extrapolate the vector meson mass, nucleon magnetic moments, form factors, charge radii and strange form factors 
HH . Finally, we note that FRR has the unique advantage that it provides a natural connection between physical results 
and quenched lattice data [4II l4ll 1^ l46j . 

In this paper, we will focus on the low moments of the GPDs. The LHPC lattice data will be used as input for 
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the extrapolation. The contribution from the disconnected quark loops has not been included in the lattice QCD 
simulation, which means that part of the sea quark contribution has been omitted [29j . While this contribution 
cancels in the isovector moments it could be very important in the isoscalar moments. Extensive investigation of 
other nucleon properties suggests that this omission might be more important at the physical quark mass than at the 
relatively heavy masses where the lattice calculations were made [32j. If this were the case, the chiral extrapolation, 
which does include contributions from disconnected diagrams, may yield a reasonable representation of the physical 
values, even for isoscalar quantities. This is indeed the case for the octet magnetic moments [U EU, for example. 
On the other hand, we are unable to quantify the error associated with this procedure for isoscalar quantities and 
in future it would clearly be preferable to be able to work with lattice simulations which include the disconnected 
contributions. 

The paper is organized as follows. In section II, we briefly introduce the chiral Lagrangian which is used for the 
calculation of the first moments. The first isovector and isoscalar moments will be derived in section III. Numerical 
results are presented in section IV and finally section V presents a summary of the results with some discussion. 

II. CHIRAL LAGRANGIAN 

The generalized parton distribution functions, H q (x,^,t) and E q (x,l;,t), are defined as 

^ e iXx (p',a'|V (~) M (t) |P ' S) = H*(x,Z t t)u(p,8)+u(p,s) 

+ E q (x, £, t)u(p' , a ) u & s ^ W 

where Q — p — p,t = Q 2 ,£ — —\n ■ Q with n the light-like vector satisfying n 2 = 0, \n ■ Q = 1. These GPDs can 
be transformed into Mellin moments/form factors by integrating with different powers of the momentum fraction, x. 
The zero-th order moments correspond to the Dirac and Pauli form factors as 



J dxx°H q (x,£,t) = F?(t), 
J dxx°E q (x,^t) = F$(t). 



(2) 



(3) 



The Dirac and Pauli as well as the electric and magnetic form factors have been discussed widely in the literature. In 
this paper, we focus on the first moments of the nucleon GPDs: 
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dxxH q {x,i,t) = A%{t) + (-2£) 2 C 2 yi), (4) 



dxxE q (x,U) = Bl, (t) - (-20 2 Cf, (*). (5) 

i 



The form factors A, B and C can be calculated with the following equation: 



i{p Wi{ti D u}q\p) = u(p 



2 B 2 0(Q 2 ) C| (Q 2 ) 
A 2,o(Q )l{fiPu} 2 ' M Q a iO- a {^Pv} + ' M Q{»Qu} 



u(p), (6) 



where the bracket, {...}, denotes the symmetrized and traceless combination of all indices in the operator. Af/v is the 
nucleon mass and p is the sum of the initial and final momenta. As one can see from the above equation, the first 
moments of the GPDs can be calculated by inserting into the nucleon states a tensor current which interacts with the 
external tensor field. 

It is convenient to define the isospin scalar and vector form factors X (X=A, B or C) with the combination of each 
quark's contribution: 

A7+ d (Q 2 ) = A1 (Q 2 ) + X 2 d (Q 2 ), (7) 
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X^ d (Q 2 )=XUQ 2 )~Xi (Q 2 ). 



(8) 



In chiral perturbation theory, the interaction between tensor current and the external tensor field as well as the 
baryon-meson interaction can be written in a series of powers of momentum of the tensor and meson fields [35] . The 
lowest order Lagrangian is 



£(0) 
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"'°uV^r 3 uhB + a s 2t0 V° u 
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where VZ V and V® are the isovector and isoscalar tensor fields. The covariant derivative D M is defined as D M = 

z5 M — D^, where £^ = = d^ L + | [u^, <9 M u] . U = u 2 is the non-linear realization of the Goldstone boson field. 
As in Ref. |35j . the parity-odd tensor interaction term is also included. In the above equation, Aa 20 is only poorly 
known and it is related to the spin-dependent analogue of the mean momentum fraction, namely < Ax > u -d- 
The 0(p r ) part of the interaction Lagrangian is expressed as 
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[D a , uW* v T 3 U + «V*,tV] (T a ^D^ + h.C. 



(10) 



where m m is defined as — iu^d^u* . The 0(p 2 ) part of the interaction can be written as 
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where the first part of this Lagrangian is the interaction between the pion fields and the external isoscalar tensor field. 

is the momentum fraction of the pion carried by quarks. Its value is less than 1 since some of the momentum 
of the pion is carried by gluons. In the preceeding equations, a 20 , a 2 ,o> ^2,oj ^2,o> c 2,o an£ l c 2,o are the low energy 
constants which can be determined from the lattice data. 

In our calculation, we also include the A intermediate baryon, which is known to be crucial in the calculation of 
spin dependent quantities. The interaction between the A and the external tensor field can be written as 
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where F 3 ^ = [D a , V 3 V ] - [Dp, V 3 V ] and F^ v = [d a ,V$„] - [dp,V£ u ]. The labels d and t in the low energy constants 
stand for decuplet and transition, respecti vely . T a are for decuplet fields which have three flavor indices (they are 
not shown explicitly, see for example, Ref. [5(| for details), defined as 



T m = A++, T U2 = A+, T 122 = A°, T 222 = A- 



(13) 



Within the framework of SU(6) symmetry, there are relationships between the octet and decuplet coefficients. For 
the isovector coefficients: 
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For the isoscalar coefficients, we have: 

s,A + + _ s s,A + 
a 2,0 — a 2,0' °2,0 

,s,A++ _ ,, s ?, S ' A+ — 

°2,0 — <J0 2,0> °2,0 — 

S,A + + s S,A^ S S,l±' s S,iA s /-| r\ 

C 2,0 — C 2,0' C 2,0 — C 2,0' i C 2,0 ~ C 2,0' > C 2,0 ~ C 2,0- l 10 7 

With the Lagrangian, we can caiculate the form factors. In the case of the lowest moments, the electric form factor 
is related to the contribution from the time component of the vector current, while the magnetic form factor is related 
to the space component contribution. For the first form factors, similar as in the case of electric and magnetic form 
factors, one can also get three form factors .42, o, B2.0 and C2,o in the heavy baryon formalism, whose relationships to 
the tensor current are 

J« ee i(p'\q 1{0 1$ 0} \p) = 3 -fA%{Q 2 ) + ^-C 2 %(Q 2 ), (16) 
J| 3 ee i(p'\q 7{ A } \p) = |^, (Q 2 ) + ||^2%(<n (17) 

j* ee l {p\q 1{0 ^ 3} \p) = ^^,o(Q 2 )(^ x Qh- (18) 

Within the heavy baryon formalism these three form factors are related to the commonly used alternative form 
factors: 

eUQ 2 ) = b 2 %(<3 3 ) + ai „(Q 2 ) - 8Ma , ( e + m„) b ''° 1q2) (20) 

ci,m = <*.<«"> + UI J 2 +MN f UQ') (2D 

III. FIRST MOMENTS OF GPDS 

In this section, we will derive the formulas for the isovector and isoscalar form factors ^2,0 , ^2,0 and C2,o- The one 
loop Feynman diagrams are shown in Fig. 1. The solid lines are for the nucleons and A. The dashed lines are for the 
7r meson and the dotted lines are for the external tensor current. Not all the diagrams have contribution to isoscalar 
or isovector form factors. From Eq. dill) , one can see that there is no corresponding interaction between the pion and 
the isovector tensor field. This is because the r matrix is traceless. For the isoscalar form factors, the contributions 
of the 7r + and ir~ loops in diagram c cancel. As a result, diagram c gives no contribution to either the isoscalar or 
isovector form factors. In Fig. 1, graphs b, d and f are of order 0(m 2 ), while e is of order O(m^). The order of 
diagram a is dependent on the moments. Its order is of 0{m%) for A20 , 0(m%) for S|.o anc l C( to tt) for C| . We 
study the isovector form factors first. For the isovector form factors, the diagram Fig. la and lc have no contribution 
and diagram le only contributes to .4^ . 



The contribution of Fig. lb, including wave function renormalisation (Fig. If), is expressed as 

. z „ So _ |l% - £Sj U k M» , ,22, 
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Fig. le only contributes to the form factor ^2,0 expressed as 

A ^ ~ 4*n*F*M N J d k W 2 ( fc) ■ ^ 

From the formulas one can see that in the heavy baryon formalism, the lowest order loop contribution to the isovector 
form factors are Q 2 independent. 

The total expression for the form factors can be written as 

/ 

Q2AQ 2 , ™l) = Z(gt + glml + (g v t + gl^l)Q 2 ) + ]T Q v 2 % (30) 

i— a 

where ^0 stands for A% , ^2,0 an d ^20- 92 o> 5tt; 9t an d 9% t are the corresponding low energy constants which are 
determined by fitting the lattice data. In particular, g% is identical to a\ , b\ and C2 , correspondingly. The other 
terms in the above equations are from the tree level Lagrangian of high order. For example, for A% oj the m 2 dependent 

term can be obtained from the interaction -tpN | [u^V^^u + uV^ v t 3 u^ (x+) 7^ U tpN, where \+ = u^x u ^ + U X <U 

and x — 2BM.. B is the chiral condensate and M is the quark mass matrix. (...) denotes the trace in flavor space. 

The Q 2 dependent term comes from the interaction ipN | [D a , [D a , u^V^ v t 3 u + uV^„r 3 ut]] 7^ D ippf. 

Since in the heavy baryon formalism, the loop contribution is Q 2 independent, the Q 2 dependence appears only in 
the analytic part. One can also fit the data versus Q 2 at fixed pion mass. In that case Eq. (|30l) becomes: 

/ 

<32,o{Q 2 )\fixed ml = Z (h\ + h^Q 2 ) + ^ £3,0 > ( 31 ) 

i— a 

where h\ and h v 2 are fitted independently for each pion mass. 

For the isoscalar form factors, Fig. la, Fig. lb and Fig. If contribute, while fig. lc and Id give no contribution. 
The reason is that in Fig. lc and Id, the 7r + and 7r~ loops cancel each other exactly. The contributions from lb and 
If are 

- z«u + m + Ss [*> 1^ .. . («» 



64tt 3 F2 J ' w 3(g) 8tt 3 F 2 J tj(k)(uj(k) + 5) 
B s — Zb s 3g ^° (^/^ 5C2bl ° [d*k Pu2{%) (33) 



192^ F 2 J w 3(jfe) 72tt 3 F 2 J u (k)(u(k) + S) 
Q. = *J „ + + £k l*k **® . (34) 

2,0 Q 4n 3 F 2 J 8tT 3 F2 7 W (fc)( W (fc)+<5) 2 

Again, the first integrals of the above three formulas are from the nucleon and the second integrals are from the A 
intermediate state. 

To evaluate Fig. la, we need to calculate the contribution of three components of the tensor current defined in 
Eqs. (p~6|) - (p~8|) . The expression for these components are 



* - -^/ ^*^(Vtf-V a '' - jSI ' <w " (t) " (t - »>* ■ <* - (35) 

^ = 3gX [ d 3 k u(k)u(k-q)[4k z (k z — q z ) — /■•■( /v - 7)] A- ■(/••- </) 



l28ir 3 F 2 J oj 2 {k)uj 2 (k-q) 
48tt 3 F2 



d A ku{k)u{k - q)[4:k z {k z - q z ) - k ■ (k - q)}k ■ (k - q)f(u), (36) 
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where /(w) is expressed as 



w(k) + uj(k -q) + S 
u}(k)[uj(k) + S]uj(k - q)[uj(k - q) + S][u)(k) + u{k - q)] ' 



(37) 



(38) 



Both the intermediate nucleon and A states are included. One can see that this diagram gives the Q 2 dependence of 
the form factors. The total isoscalar form factors , £>| and C| can be written in the same way as Eq. (|30|) . The 
corresponding low energy constants can also be determined by the lattice data. 



sf 0.4 




FIG. 2: The form factor A^X d versus Q 2 at each pion mass and versus pion mass at Q 2 = 0. The dashed and solid lines 
correspond to the global and separate fit to the lattice data, respectively. 



IV. NUMERICAL RESULTS 



In the numerical calculations, the parameters are chosen to be D = 0.76 and F = 0.50 (gA = D + F = 1.26). 
The coupling constant C is chosen to be — 2D, as estimated by SU(6) relations — which gives a similar value to that 
obtained from the hadronic decay width of the A. Here the finite-range regulator is chosen to take the dipole form 

u(k) = , (39) 

(l + fc 2 /A 2 ) 2 

with A = 0.8 GeV. The regulator has been applied in our previous work on nucleon mass, magnetic moments, form 
factors, charge radii, etc. The other two parameters in the Lagrangian and AaJJ & r e chosen to be 0.7 [5l| and 
0.21 [35j, respectively. 

All the lattice data had been transformed to a scale of /i 2 = 4 GeV 2 . We first study the isoscalar form factors. 
The form factor A%~Q d is shown in Fig. 2. The lattice data are from Ref. [29|. The solid lines correspond to fitting 
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FIG. 3: The form factor £>2 Q d versus Q 2 at each pion mass and versus pion mass at Q 2 = 0. The dashed and solid lines 
correspond to the global and separate fits to the lattice data, respectively. 

the lattice data with Eq. (|3"TT) at each pion mass separately. The dashed lines correspond to fitting the lattice data 
with Eq. (|30|) for all the pion masses. We can see from the figure that the difference between the solid and dashed 
lines is relatively small. Both of them can extrapolate the lattice data very well. At each pion mass, the lattice data 
shows a dependence on Q 2 which is almost linear. In fact, the result of the Fig. 1 is Q 2 independent except for Fig. la, 
which gives a little curvature to the line. It is interesting that the momentum dependence of the first form factor 
is not like the electromagnetic form factor which has a dipole behavior. For the dependence, it is not like the 
electromagnetic moment either. For example, For the magnetic moment, Fig. la gives the leading order of 0(m n ) 
contribution. For , the leading order is of 0(m 2 ). This is because in this case, Fig. la are of 0(m^). Fig. lb 
and Fig. If are of 0{m 2 ! ) which is the same for the magnetic form factor. Therefore, the magnetic form factor has a 
large curvature at small pion mass. At Q 2 — 0, A^\ d increases with increasing pion mass. 

The form factor B\^ d is shown in Fig. 3. The solid and dashed lines have the same meaning as in Fig. 2. Again, 
the lattice data can be described very well. Similarly to A% g d , has little curvature with increasing Q 2 . The 

form factor is not sensitive to the pion mass. For this form factor, all the diagrams in Fig. 1 are of 0(m 2 ) including 
Fig. la. At the physical pion mass, the value of B 2 ^ d at Q 2 = is 0.497 ± 0.089. 



From Ji's sum rule [52| , the relationship between the moments and the quark contribution to the total nucleon spin 

is 

Ju+d = \[A s 2fi {Q 2 = 0) + Bl (Q 2 = 0)] = \B s 2fi {Q 2 = 0). (40) 

With the extrapolated value of £?| (Q 2 = 0) at the physical pion mass, we find J u +d = 0.249 ± 0.045. Since the total 
spin of the nucleon is 1/2, it is interesting that only 50% of the nucleon spin is carried by quarks. This is consistent 
with studies of the evolution of the total quark angular momentum from a scale typical of a valence-dominated quark 
model [H . 

In Fig. 4 we show the form factor C^"^. Once again the Q 2 dependence shows little curvature. However, this time, 
the m-x dependence has a visible curvature at the physical pion mass and Q 2 — 0. This is because, unlike A 2 ~\f and 
^2 t •> Fig- la has a leading nonanalytic term of order 0(m n ) for C^ 1 ^- The absolute value of decreases with 
increasing pion mass at Q 2 =0. 
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FIG. 4: The form factor versus Q 2 at each pion mass and versus pion mass at Q 2 = 0. The dashed and solid lines 

correspond to the global and separate fits to the lattice data, respectively. 
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FIG. 5: The form factor A% d versus Q 2 at each pion mass and versus pion mass at Q 2 — 0. The dashed and solid lines 
correspond to the global and separate fits to the lattice data, respectively. 



Having discussed the isoscalar moments which describe the total contribution of the u and d quark, we now turn 
to the isovector part. The isovector moments describe the difference of the u and d quark contributions. In this 
case, diagram a and c in Fig. 1 have no contribution and, as a result, the isovector form factors of the loop are Q 2 
independent. The form factor A^q 1 is shown in Fig. 5. Just as for the isoscalar form factors, the diagrams of Fig. lb, 
Id and If are of 0(m 2 r ). Because Fig. le is of 0{m%), the result is not sensitive to the choice of Aa^ - For example, 
at Q 2 = 0, the values of A% are 0.156 and 0.153 for Aa^Q = 0.21 (phenomenological value) and 0.144 (fit value), 
respectively [H, |3(|. For the Q 2 dependence, from the figure one can see that the lines are well described by a linear 
dependence on Q 2 , which comes from the choice of tree level contribution. 
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FIG. 6: The form factor B%^ d versus Q 2 at each pion mass and versus pion mass at Q 2 = 0. The dashed and solid lines 
correspond to the global and separate fits to the lattice data, respectively. 
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FIG. 7: The form factor C% d versus Q 2 at each pion mass and versus pion mass at Q 2 = 0. The dashed and solid lines 
correspond to the global and separate fits to the lattice data, respectively. 



In Fig. 6 we show the isovector form factor £>2 d - At Q 2 = 0, it increases with increasing m 2 . At the physical 

pion mass its value is about 0.433 ± 0.071. B^ d is close to which means that the u quark is dominant for the 

proton spin while the d quark gives little contribution. The form factor C\ is shown in Fig. 7. The value of C\ is 
around zero and not sensitive to either Q 2 or m 2 . 



11 



TABLE I: Low energy constants and moments at physical pion mass. The results in the table are for the linear fit where the 
momentum dependence of the tree level term is up to Q 2 . For the dipole fit, the first moments at Q 2 = are about 10% — 20% 
larger. The results of LHPC are also listed in the last column. 





ff2,0 


^(GeV" 2 ) 


St(GeV- 2 ) 


g.,i(GeV- 4 ) 


02,0 (0) 


5 2 T C (0)[29] 


AU + d 

"^2,0 


0.489 


0.152 


-0.266 


0.019 


0.440 ± 0.033 


0.520 ± 0.024 


fju+d 
"2,0 


0.496 


0.005 


-0.339 


0.268 


0.497 ± 0.089 


0.425 ± 0.086 


rtu-\-d 
C 2,0 


-0.152 


0.102 


0.119 


-0.159 


-0.217 ±0.103 


-0.267 ±0.062 


AU — d 

•^2,0 


0.207 


0.010 


-0.139 


-0.009 


0.156 ±0.020 


0.157 ±0.010 


ryU — d 

°2,0 


0.526 


0.158 


-0.298 


0.111 


0.433 ± 0.071 


0.430 ± 0.063 


riu — d 
°2,0 


-0.001 


-0.009 


-0.025 


0.062 


-0.001 ±0.050 


-0.017 ±0.041 



The values of ^2,0 , &2,o and £2,0 at the physical pion mass and Q 2 = are shown in Table I. Using Eqs. ([T^|) - ([2T|) . 
one can extract the form factors ^2,0, -62,0 and C^o- In particular, at Q 2 — 0, ^2,0 = ^-2 .0 , £>2,o = ^2,0 ± -62,0 and 
£2,0 = ^2.0- Compared with those in Ref. (29j extrapolated with the formulas in Ref. [35| . our numerical results are 
close to their results which used dimensional regularization. With the values for the isoscalar and isovector moments 
in the Table, one can easily find the u and d quark moments. For example, for the moment A2, 0, the contribution 
from the u quark is about as twice large as that from the d quark. This is because in the proton, there are two u 
quarks and one d quark. The values of the quark spin are J u — 0.233 ± 0.04 and Jd = 0.016 ± 0.04 . The u quark 
dominance for the moment 62,0 can be understood from the naive quark model for proton, where the u quark spin is 
as four times large as the d quark spin. For C2,0j it seems that the u and d quarks yield almost the same contribution. 

0.3 -p — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 



0.2 



0.1 




0.0- 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Q 2 ( GeV 2 ) 

FIG. 8: The form factor versus Q 2 at m n = 758 MeV. The solid and dashed lines are for the linear fit as in Fig. 5. The 

dotted line is for the dipole fit. 

We should mention that our extrapolation results depend on the choice of the tree-level contribution. In the previous 
fit, the order of the momentum dependence of the tree- level contribution is up to Q 2 . This is supposed to be valid at 
low momentum transfer. We know that for the electromagnetic form factors, the momentum dependence has a dipole 
behavior with mass parameter around 0.71 GeV 2 . In the case of the axial form factor, which is perhaps more relevant 
for the spin distribution, the phenomenological form factor is a 1 GeV dipole. Since the actual lattice data extends 
over such a broad range of Q 2 , a pure linear dependence on Q 2 is difficult to justify. Therefore, we have carried out 
another fit with a modified dipole form. The Q 2 dependence of the first moments, i.e. Eq. (f3"Tj). is now changed to 
the following expression: 



SloiQ 2 ) 



fixed m\ 



Z{h\ + h v 2 Q 2 ) 
(1±Q 2 /A 2 ) 2 



T, r " 



(41) 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Q 2 ( GeV 2 ) 

FIG. 9: The form factor g d versus Q 2 at — 758 MeV. The solid and dashed lines are for the linear fit as in Fig. 6. The 
dotted line is for the dipole fit. 

where A is chosen to be 1 GeV. While a pure dipole form would be prefered phenomenologically, we note that the 
data from Ref. [29j tends to give much harder form factors for the axial channel than those found in nature. The 
reason for this is not understood but our modified form allows room to fit the lattice data while still including some 
physically reasonable Q 2 dependence. 

We show the momentum dependence of A% Q d and B^ g d at m, = 758 MeV in Fig. 8 and Fig. 9. As a comparison, 
the solid and dashed lines with the previous linear fit are also shown. The dotted lines correspond to the fit with the 
modified dipole form at tree level. The lattice data can be reasonably described with both forms. For £>2 g d , compared 
with the linear fit, the dipole fit gives a larger value at Q 2 = 0. The situation is similar for the other pion masses. 
The curve of mass dependence of B\ in Fig. 6 is shifted up in the dipole fit. As a result, the moment at the physical 

mass is increased from 0.433 to 0.53. For A% g d , the difference between two fits is not as large as for B% ~^ d - This is 

because the lattice data for A^~ d at Q 2 — impose a strong constraint in that case. At the physical pion mass, its 
value changes from 0.156 to 0.17. 

The dipole fit makes the absolute value of all the first moments larger at Q 2 = 0. For .4.2,0 > the difference is about 
10% with the help of the lattice data at zero momentum. For £12,0 and C2 o d > the difference is about 20%. The value 
of C2 Q d is still around zero. The difference between these two separate fitting procedures provides some indication of 
the systematic error in extracting this important physical information from the lattice data. 

V. SUMMARY 

Chiral perturbation theory with finite-range-regularization has been applied to the problem of extrapolating lattice 
QCD data for GPD moments to the physical pion mass and zero momentum transfer. For the isovector form factors, 
the one loop contribution is of 0(m 2 ). For the isoscalar form factors, A% g d and £> 2 Q d ; the leading order is of 0(m 2 ), 

while for C% g d , the leading order is of ©(m^). The lattice data were fitted both globally and separately at each pion 
mass. At Q 2 = 0, the dependence of the first moments (except C 2 Q d ) does not show a big degree of curvature at 
small pion mass, which is quite different from the zero-th moments (electromagnetic form factors). Overall the level 
of agreement between the extrapolated results obtained using dimensional regularization and FRR is satisfactory. 

The Q 2 dependence, especially of B2.0, is mainly determined by the tree level Q 2 behavior. In our first calculation, 
we retained only the tree level terms up to Q 2 , that is a linear dependence. On the other hand, the data has been 
determined over a wide range of Q 2 , up to 1.2 GeV 2 . Phenomenlogically one knows that the physical form factors 
exhibit a considerable variation with Q 2 over such a range - usually described by a dipole form. In order to test the 
sensitivity to this problem, we have re-done the fits to the lattice data using a modified dipole form. In this case the 
first moments at Q 2 — are about 10% — 20% larger than those in the linear fit. We regard this as a measure of one 
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of the main systematic errors associated with extracting information about the quark angular momentum. 

In our calculation, the A contributions have been included explicitly, because of its potential importance for spin 
dependent quantities. For the nucleon spin, the inclusion of the A changes the result by 10% [7|. However, for the 
extrapolation of the first moments of GPDs, if the A is not included the low energy constants are different. Since 
they are adjusted by fitting the lattice data, the extrapolated moments at the physical pion mass change by less than 
about 5% when the A is included. 

Using the extrapolated moment 62.0, we can extract information concerning the quar k contribution to the nucleon 
spin. Our results are consistent with the current JLab and HERMES experiments [53|,l5J|. On the other hand, the 
experimental errors there are quite large at the present time. The comparison with model predictions, such as those of 
Refs. 0,111 , is also quite satisfactory, especially given the systematic errors associated with the Q 2 dependence of B 2 ,o 
(discussed above) and the dependence of the proton axial charge on lattice volume. As for the other two moments, 
the u quark contribution is twice as large as that for the d quark in ^2,0, while they give similar contributions to 
C2,o- In comparing with data, it is important to note the possible sensitivity of isoscalar quantities to disconnected 
graphs which were not included in the LHPC simulation. In the next few years we can look forward not only to more 
accurate data on GPDs but lattice QCD data on larger volumes and at lower quark masses and, ideally, including 
disconnected terms. All of this offers enormous promise towards unravelling the proton spin problem. 
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